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Abstract 


This document is a final technical report on a National 
Aeronautics and Space Administration research grant titled 
"Three Body Recombination in Supersonic Flow", NSG 1314, June, 
1976 through July, 1978, from NASA Langley Research Center. 

The investigators proposed to look for diagnostic techniques 
and experimental evidence to support new approaches to the 
turbulence problem using the shock tube and to undertake 
related theoretical studies as well. Evidence to support a 
kinetic theory of turbulence has been found. An overview of 
these and other achievements to date is given herein. 
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Table Captions 


Table I 


The Performance of Pressure Loaded Arc Driven 54 . 
Shock Tubes. See text for the explanation of 
symbols. 
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Figure Captions 


Figure 1 


Figure 2 


The- growth of the turbulent intensity of bursts 13 

in an unstable detonation front. The Reynolds 

number is based on the displacement- thickness as 

measured from an interferogram of the front from 

2 1/2 

reference 9. The quantity (AS ) is the root 
mean squared averaged value of the fluctuations 
about an' average fringe shift in the turbulent 
burst. The quantity I is defined in the text. 

The dashed straight line is a fit to the four 
data points indicated which shows a slope 
(d(log I)/dRe = 3.2 x 10” 3 ) consistent with 
theoretical expectations. 

Overall Shock Tube Flow. This is an (x,t) 16 

diagram. A partition separates ' two regions in 
the shock tube with high pressure gas occupying 
x < o at t < o and the gas originally at low 
pressure occupying x > o at t < o. Examples 
of the various discontinuities and wave phenomena 
in the shock tube subsequent (t>o) to the rupture 
of the partition at x = o, t = o are shown. 

(S) is the shock wave; (C) is the contact surface; 

(T) is the tail of the expansion fan (in this 

case (v^-a^) > 0); (H) is the head of the expan- 

sion fan. (See text.) 
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Figure 3 


19 


Figure 4 


Figure 5 


The Test Section. Some of the pressure gauges 
are indicated along the top of the tube. The 
optical path through the test section is shown 
in its Michel son-Mor ley interferometer mode. A 
single fringe from the interferometer is shared 
by two light pipes which provide, with associated 
photo-electronics, signals FI and F2 respectively. 

The insert shows a sample of calibration data for 
F = F1-F2 confirming fringe shift cancellation for 
unsteady irregular flow where the correlation 
length scale is large compared to the separation 
between the light pipes. 

Calibration: Shock Front Pressure Ratio. Here 22 

AP = {(V P l>msas ' <V P lW ! ( V P l> C al is 

the pressure ratio calculated from the measured 
Mach number of the shock wave (S in Figure 1), 

M , and (P 0 /P 1 ) is the pressure ratio across 

the shock wave as measured in the test section. 

Calibration: Contact Surface Time of Flight. 24 

T M is the measured time interval between the 
c ,M 

arrival of the shock wave at x = +76.4 cm in the 
test section and the appearance of rapid irregular 
density fluctuations {i.e., at Fl) as indicated 
by one of the light pipes in Figure 2. This is 
schematically indicated in the sketch where P 
shows pressure versus time and the arrival of 

(vi) 



Figure 6 


Figure 7 


the shock wave and D shows the output of the 

interferometer with the amplitude modulation 

resulting from fringe shifts. A sample of the 

data, with phototube output (arbitrary units) 

versus time (50 usec/div.), is displayed as an 

insert. T is the calculated value of the 
c, c 

interval, using the shock wave's measured Mach 
number, M , and assuming ideal gas flow, between 
the arrival of the shock wave and the arrival of 
the contact surface at the same location. 

Calibration: Lambert- Beer ' s Law Absorption for 25 

O 

NO 2 . Relative absorption of blue light (4350 A) 
for several concentrations of NC> 2 appropriate to 
our operating conditions is shown. The straight 
line is a best fit to (E/E q ) = 10 a where 
a = -(0.208 ± 0.004) [N0 2 ] . 

Examples of Large Fluctuations in Pressure. In 27 

each photograph, the scales are: vertically, 

pressure (.14 atm/div) ; horizontally, time 
(.5 msec/div) . The middle trace (ignore the 
uppermost trace) is pressure at x = +76.4 cm 
and the bottom trace is pressure at x = +122.1 cm, 
further away from the partition. Each photograph 
shows one or more large fluctuations at x = +76.4 cm 
subsequently diminishing in deviation at 
x = +122.1 cm. (a) M = 2.2, (ReVcm) = 3.9 x 10 5 ; 
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Figure 8 


28 


Figure 9 


(b) M = 2.2, (Re^cm) = 4.4 x 10^; (c) M = 2.0, 

S 3 

(ReVcm) = 2.7 x 10 5 ; (d) M « 2.0, (reVcm) = 

b 

3.8 x 10 5 . 

An Example of Total Gas. Density Sudden Coherence 
Associated with a Large Pressure Fluctuation. The 
top trace is the pressure gauge's output. The 
bottom trace is the (simultaneous) output of 
F (= F1-F2) from the interferometer located at 
the same position in the test section as the 
pressure gauge. The traces shown are as displayed 
on an X-Y plotter which records the output of the 
Biomation 8100 ADC. The arrow indicates the 
large fluctuation in pressure. 

An Example of Reactant Gas Density Distortion 29 

Associated with a Large Pressure Fluctuation. ' The 
top trace is- the pressure gauge's output. The 
bottom trace is the (simultaneous) output from the 
phototube which monitors the absorption of light 

O 

at 4350 A along the optical path through the test 
section. The decreased N0 2 density shows up as 
decreased absorption and increased negative signal 
E at the phototube. The traces shown are as dis- 
played on a computer reconstruction of the output 
of our analog to digital converter. The arrows 
indicate the coincident fluctuations. 
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Figure 10 


30 


An Example of a Structured Correlation Function 
during a Large Pressure Fluctuation. A segment by 
segment quadratic fit to the pressure history p (t) 
at the test section's window is made after the data 
are divided into overlapping 200 ysec segments. 

Segment S is centered on the pressure fluctuation. 

The instantaneous value of the pressure during 
a segment, p(t), is compared with the value obtained 
from the fit for that segment, (t) ; p' (t) is 
given by p* (t) = p(t)-p f (t). The separation 
between the data shown above. A, is 200 usee. 

R{t) is defined in the text. 

Figure 11 Fluctuation Dynamics. (a) Fluctuation Width vs 32 

Time Lag. The duration of the fluctuation as 
evidenced in the pressure history is given by twp. 

The simultaneous duration of the fluctuation as 
evidenced in the interferometer data (F) is given 
by twi. The time lag between the arrival of the 
contact surface and the appearance of the fluctu- 
ation at the optical window in the test section is . 
tD. (b) The Inverse of the Time Lag vs Pressure 
Deviation Strength. The peak value of the pressure 
in the fluctuation is p_. The average value of 
the pressure immediately before and after the 
fluctuation is p^. The straight line is a least 
square fit (correlation coefficient = 0.997)' to 
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w = 1661 Z 


-1 -3 -1 

- 1405 in which w = tD x 10 sec 
and Z = Pr,/p T • The data in (b) are taken from 

x3 Xj 

the 3 2 fluctuation (see text) . 

Figure 12 Fluctuation Dynamics: Velocity. With the velocity 33 

of the contact surface, v^ (computed from the 
measured value of M ) , and the measured velocity 
of the pressure fluctuation as tracked' through two 
consecutive pressure gauge histories, v^, then 
Av H v 3 ~v l . 

Figure 13 Fluctuation Dynamics: Reynolds Number. (See text.) 3M 

The time lag between the appearance of the contact 

surface and the appearance of the large pressure 

fluctuation at the test section's window is tD. 

The Reynolds number per unit centimeter is computed 

from the measured Mach number M for conditions in 

s 

the driver gas at the contact surface. 

Figure 14 Boundary Layer Model: log (1/tD) vs Re. The 38 

Reynolds number is derived from a boundary layer 

model for the contact surface; the quantity tD 

is as defined previously. The dashed lines are 

least squares fits with correlation coefficients 

r D = 0.992 and r D = 0.998. The average slope 
B 1 3 2 

for the two lines is m Av = 2.7 ± 0.7 where 
m = d{log (1/tD) }/d (re) . 
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Figure 15 


Shock Heated Driven Gas Flow Velocity vs. 
Driver /Driven Gas Density Ratio. Solid line: 


45 


Y-j_ - 7/5, y 4 - 5/3. • Dashed line: = y 4 = 5/3. 

Figure 16 Shock Wave Speed (Arbitrary Units) vs. Driver 47 

Energy Density. (See text for details.) 

3 4 

o, Menard ; o, Gruszczynski and Warren ; 

4 

A, Camm and Rose ; x, present measurements. 

Figure 17 Mach Number vs. Initial Driver Tube Pressure. 50 

A, P.^ = 0.1 Torr; +, p^ = 0.5 Torr; o, p.^ = 1.0 Torr; 
o, p.^ = 2.0 Torr; x, p 1 = 0.4 Torr. All data shown 
are for a He/Ar set-up. 

Figure 18 .Electrical Conversion Efficiency vs. Initial Driver 51 

Tube Pressure. The data shown in Figure 3 have 

been averaged over p.. . 

1 > 

Figure 19 Absorption Coefficient vs. Temperature at a Fixed 62’ 

Gas Density. (1) n « 10^ cm (2) n — 10"*"^ cm 

ti\ -I r\l-7 — 3 f a \ . « 1 6 —3 

(3) n == 10 cm ; (4) n = 10 cm 

Figure 20 Absorption Coefficient vs. Temperature at a Fixed '64 

19-3 

Electron Density. (1) n = 10 cm ; (2) n = 

e • e 

rt 18 —3 i nl7 —3 , n 16 —3 

10 cm- . ; (3) n = 10 cm ; (4) n = 10 cm 

e e 

Figure 21 Absorption Coefficient vs. Electron Density at a 65 

Fixed Temperature. (1) T = 12,000°K; (2) T = 18,000°K; 
(3) T - 24 , 000 °K. 
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I. Introduction 

We have been motivated to investigate alternatives to 
the phenomenological approach to turbulence by the well- 
known deficiencies in that approach. These deficiencies 
include: (1) the use of empirical constants not derivable 

from the fundamental constants of nature; (2) the inability 
to consistently explain the existence and systematics of large 
coherent structures in "steady" turbulence; (3) the inadequate 
characterizations of transport effects observed in reacting 
transitional flow, especially, flow with exothermic processes; 

(4) the incomplete statistical description of the randomly 
fluctuating gasdynamical variables. Generally stated, we have 
sought to create and understand turbulent environments where 
phenomena might be observed which would allow us to distinguish 
and evaluate new approaches to turbulence. In parallel with 
these efforts, we have sought to extend the new approaches 
into regimes which are accessible to our measurement capabilities 
This report reviews the principal specific results within that 
general context from our beginning two year studies of turbu- 
lence with NASA support. 
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II. Reaction Rate Distortion in Turbulent Flow 

Consider the irreversible chemical reaction A + B P 
in a turbulent gas from a classical point of view. Impact 
parameters between particles 5. and k are b^, the separation 
between their centers, and e, the angle between their original 
and final velocities in the center of mass. Each particle 
has a characteristic mass m and a characteristic internal 
energy E^. Let's assume that: (1) interactions are unimpor- 

tant; (3) only binary collisions need be considered. Mass 
and total momentum are to be strictly conserved in all 
collisions. Energy is also conserved; however, the product P 
is produced in an excited state (with excess internal energy 
E p *) from which it very promptly relaxes to its characteristic 
internal state through some unspecified noncollisional 
(radiative) process. We shall apply the approach used by Tsuge^" 
to this problem as we look for an explicit influence of 
turbulence on the apparent reaction rate. 

2 

Using the -Klimontovich N-particle distribution function 

O 

with the six dimensional phase space coordinate Z, f (Z,t) s 
N (s) 

I <$[Z-Z V ' (t)J, the Master Boltzmann Equation for A may be 
s=l 

written as: 



II- (1) 
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where & = A, B, then P and the function "F" is the fraction 
of collisions between A and B which result in formation of 
P. The equation for B is the same as II- (1) with the B label 
replacing A: 


, ‘3 . r 3 ® 

3t Vfi 3x r B 


(B-hA above) }- 


b AB db AB dEd3v A 


II- (2 ) 


For the product P: 




O 




II-(3) 


Equation II- (1) simply says that the rate of change in the 

3 3 

number of molecules in a phase-space volume element d xd v, is 
equal to the sum of the differences of the rate of molecule A 
into and out of the specified volume element due to collisions 
between: A and A; A and B; A and P; less the rate resulting from 

the formation of P. Equation II- (2) has the same meaning for 
molecule B. Equation II- (3) is different since there are no 
collisions involving P which can cause the reaction; however, 
a fraction of the collisions of A and B do produce P. The 
inclusion of a delta function is to insure that only those 
collisions which conserve momentum are allowed and to guarantee 



Jo M 


4 . 


that the resulting new P molecules will be in the specified 

3 3 

volume element d xd v around x and v^. 

P 1 

O 0 ^ A 

Following Tsuge, we can define: f = f; f^f B = ^ab^ Z A ,Z B^ 


f s i s = s s[V z a <Sl (t)1 V z a <S> (t)1 =_y [V z a ,(t > ] 6 ! V z 

S— 1 


(s] 


S=1 


St^S 


A 


N 


+ 5 ( z a- z fi> z , 6[z a" z a <S> lt) 1 ” f AA (z a' z a> + « (z a"-V f A <z) • 


S=1 


Turbulent-like effects are obtained without ad hoc assumptions 

from situations where, defining ip(Z,A) 5 f : 11 (Z, Z) -f (Z) f (Z) , one 

finds a history with \fj(Z,Z) ^ 0. The importance' of this approach 

3 

to turbulent , flow without reactions has already been discussed. 

By taking the averages of equations II- (1), II- (2), and II- (3) and. 
then, in turn, multiplying each equation by a .moment function 
a(v £ ) = and integrating over d 3 v £ , we can add the three 
resulting equations and find (using standard notation) : 


l 3t 


j3 9 

w i + 

9x 




S&n o 

( -TE + ^ f V i 


— j; /l. 
k,JT2j 


(a i + “k~ < V a k ) . [ f V z k> la^a^a 3 ^} 


+ 1 a p (F) [ V B +<I> M <Z a , Z B ) ] « <» a V m B V I "pV ^ VSs^V^AB 


" <“ fl + “B> < F) l a W^AB 


(t) 1 


II- (4) 
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where k also stands for A, B, or P, and d 3 K,„ « I v„ -v„ I b, _db,, _ 

AB 1 A B 1 AB AB 

Now we can define, in the usual manner and again using 

3 

standard notation, the hydrodynamical quantities of interest: 


de. 


n £ 


V 3 v i' p m £ n Jl' u «. 


5 l £ i a3v S.' S ° | £ “t m tV 


h = u sT u = n 


l, 1 . 


,3 1 f-*- 

(v 0 -u)f £ d v £ = — c £ 

V. 


f s. a3 . v i' S « 


v £ -u; 


n « Z n £ , E = (E n p E„)/n, 3KT/2 = {Z m 


Z _ 

P. . = m 0 
13 Z 


A JJ A J 


A 


■ l £ 


C £ 2 f £ d 3 v £ )/2n; 


4 c h a3v r p ij = £ p ij ;ana q i = ? j 


z 


A m.Cf e 

'A 1 + V f » dS f 


When a. = m„ and ap = m p v^ Equation II- (4) becomes, using the 


Z Z Z Z Z 

definitions above: 


+ _?_ ( Pu r ) « 0 and 3 ( gg- ) - + [P ir + pu x u r ]' = 0 
3 t „„r VK dt 3x r ir 


ax 


respectively since mass and momentum are strictly conserved. 

Thus, these conservation equations are unchanged from their 

usual appearance even for inelastic collisions. When 
1 2 

a £ ” "2 m £ v £ + E £ r suiranat i° n term on the right hand side of 

equation (4) is zero due to strict conservation of energy in 
the elastic collisions. The remaining terms on the right hand 
side are: 
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f ,-{m v +m v ) 2 - . ^ 

— + Ej (F) [f A f B +$ AB (Z A ,Z B ) ]d\d J v n d J K 2 


2m r 


A B AB 


m A V A +I VV A A o ? 

( ^ + E a +E b ) (F) [f a f^ flR (Z a ,?. R )]dXd\d K ! 


"A B t AB A' B j 


B AB 


I 


tWV (F) 'WWVV ^VVXb 


where \i = m A m B / O^trag) and ^ = ! ^a“^B * * Evaluating the left hand 
side of the equation (4), we achieve 



[~ nkT + j pu 2 +nE] + 

Sx 

p :b a-V? v2] (f) 1W 


• [|> pu 2 u r +| nkTu r +P ir u 1 +q r +nEu r ] 


- i _ 2 

Of course, the reaction is prohibited when Ep-r-E^-Eg >' j V* 

If E p > E a +E b but E p =E A -E B < J uV* 2 , then part of | jiV 2 is used 
to increase the internal energy of the system and the remainder 
is released in the form of heat. If E p < E A +B B , then heat is 

i 2 

produced in the amount (E A +E B + j vT* -Ep) with each reaction. 
Equation II- (5) is consistent with these expectations since the 
energy dependence is' contained in !, F". 

Our new result is contained in the term with a contribution 

A 

from ^ AB 2 g) • Equation II- (5) asserts that, when the reaction 
is allowed (F ^ 0), it can be affected by turbulence in a direct 
manner which is distinguishable from that arising from ordinary 
turbulent diffusion. A new apparent reaction rate fraction, 

F 

• App 


, can be defined: 
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(F) [f A f B + ^AB^ Z A ,Z B } ] and F App F = ^AB^A' a 

A ? A 


f A f B 


f f 
A B 


Thus, there is an explicit influence of turbulence on the rate 
of production of P? the reaction is distorted in all cases 

A 

where the condition (^ AB /f A fg << D is violated. 


n > 
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III. Turbulent Bursts in a Shock Tube 

Turbulent discontinuities in shock tubes are usually 

thought to be boundary layer-like insofar as transition to 

4 

turbulence is concerned. However, theoretical treatments 

which rely on the techniques from linear stability formulations 

have seemed to be inadequate for the analysis of instabilities. 

5 6 7 

in contact surfaces and detonation fronts. Tsuge’s treatment, 

by contrast, is not restricted to the small perturbation regime; 

macroscopic turbulent effects are shown to evolve from a kinetic 

theory with reduced microscopic (molecular) chaos. Furthermore, 

transition related turbulent bursts in a Blasius boundary layer 

7 

experiment are successfully interpreted. Here in this paper, 
we suggest an extension of Tsuge's approach to a wider class 
of unstable boundary layers, viz., to analogous unstable 
gasdynamic discontinuities in shock tubes. 

When one focuses on those fluctuations which experience 
maximum amplification and identifies these fluctuations as the 
experimentally observed fluctuations, the theoretical solutions 
dictate the familiar form 

I = exp t<j> (Re, A) ] . III-(l) 

7,8,9 

Solutions of equation III- (2) are well known. The Reynolds 

number has its usual meaning. The expansion parameter X has the 
units of frequency. The quantity I is an explicit function of 
turbulent intensity, I- 5 Q/Q q , where Q is the observed turbulent 
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intensity and Q q is a minimum intensity level which must be 

exceeded in order for fluctuations to grow- The function <j> 

has no empirical constants; all constants are derivable from 

boundary conditions. Therefore, this formulation is non- 

10,11 

phenomenological and has unique physical consequences. 

If one considers the solutions of equation II- (1) in 

the form of a family of curves of I' vs Re (where I* = log I) 

at selected values of the nondimensional ratio Xv/u^ (v is the 

kinematic viscosity and u^ is a characteristic free stream 

boundary layer velocity) , systematic behaviors emerge. The 

2 2 

curves show: a first regime with (dl'/dRejX), (d I'/clRe }>0; 

a middle regime which is essentially linear, with (dl’/dRe) >0, 

(d^i'/dRe 2 ) = 0; and a third regime with (dI’/3Re)>0 initially 
2 2 

and . (d I’/dlte )<0. The linear middle regime has a slope which 

we shall name m and an intercept on the Re axis which we name 

Re n - Extrapolated to I' = 0, the middle regime gives, at 

Re = Re n , Q = Q . The third regime has a maximum in I*. of 

value I' at a Reynolds number which we now name Re . 

mx p 

The extension of this kinetic theory of turbulence to 
any given set of experimental data on unstable boundary layer 
flow with turbulent bursts is therefore straightforward. If 
there is a linear regime in a plot of log Q ( exp ) vs Re ' then 
the slope of the line should be well estimated by a m appropriate 
for the Reynolds number range being treated. If there is a 


maximum in the Q* ( eX p) (= log Q^ eX pj vs Re trajectory, then 

the value of Re at Q' , , . „ is a measure of Re . One obtains 

(exp) ,mx p 

2 

from the measure of Re an estimate of X (from Xv/u,), which 

P a 
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can be called the characteristic frequency of the instability, 
and in turn an estimate of Re n * Since in the linear regime - 
the experimental data will fit log Q = b + (m) (Re) , the estimate 
of Re^ gives an estimate of Q q by 

log Q = log Q = b + m * Re . Ill- (2) 

Re=Re ° n 

n 


For example, for Shubauer's experiment (as reported in 
reference 1) the solutions to equation III- (1) in the appropriate 
Reynolds number range show: a slope m in the second regime 

independent of X, i.e., m ^ m(X) and m E (dl'/dRe) = 

_3 

(3.310.2) x 10 ? a systematic dependence of the characteristic 


wavelength on Re from the third regime which is given by 

(Xv/uij) = (4.6±0.3) x 10" 5 x exp[-(0. 34±0.02) x 10~ 3 Re ]; a 

P 

correlation between Re and Re which is found to be 

n p 

Re = (0.2910.06) + (0 . 44±0 . 02) Re . In the data with turbulent 

n p 

bursts, a linear regime is found with slope m ( eX p) = <3Q'/dRe = 

_3 

3.5 x 10 and a Q'-axis intercept b( exp ) = -11.45. A maximum 

in the Q- 1 vs Re plane occurs at Re = 3.19 x 10 3 ; this gives 

P 

2 —5 

(Xv/u^) = 1.55 x 10 from the expression above in agreement with 

the results in reference 1. In addition, using the Q' axis- 

intercept and the relationship between Re and Re (which provides 

P 

3 

= 1.4 x 10 ) given above, our estimate of Q q from equation 

_ ' _3 

III- (2) is Q = 1.5 x 10 . This corresponds to Q /Q = 120, 

o r mx o 

i.e., a required minimum turbulent intensity level which is 120 
times smaller than the observed maximum intensity in the turbu- 


lent bursts. 
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Consider next a curved unstable detonation wave moving 
in a shock tube with constant velocity u q . The non-one 
dimensional nature of this discontinuity is assumed (as is 
now conventional) to have resulted from the shear induced by 
the curvature related dynamics. Let's model the turbulent 
front by an apparent boundary layer which is transverse to and 
symmetric about the axis of propagation. That is, its smallest 
dimension is at the center of the tube? its largest dimensions 
(small with respect to the radius of curvature of the front) are 
at the shock tube's walls. This apparent boundary layer is. 
perpendicular to the "normal "boundary layer at the shock tube's 
walls. The local Reynolds number in it can be defined using a 
measured local characteristic length, an overall characteristic 
kinematic viscosity, and an overall apparent characteristic 
"free stream" transverse boundary layer velocity u^ » u q (a /8) 
arising from the curvature where a is the total angle subtended 
by the curved front. 

12 

White ' s turbulent detonation wave experiments are amenable 
to this extension. Measurements of turbulent fluctuation levels 
and local characteristic lengths are obtained from fringe shift 

13 

fluctuations in his interferograms using well-known procedures. 

Specifically, for a self-sustaining detonation into 2H 2 + 0 2 + 2C0 

with u q = 2.19 km/s and a = 0.047 rad, turbulent bursts are 

observed for which m, . * 3.2 x 10” 2 3 and Re = 208. This 

(exp) p 

2 —5 

provides (Xv/u d ) = 4.3 x 10 . . In addition, the fit to the 
linear regime and the calculated re n (= 92) give a ratio of 
maximum observed intensity to minimum required intensity 
Qj^/Qq = 2. Figure 1 shows these results. 
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Thus might a wide variety of discontinuities with- 

similar nonlinear instabilities be treated. We notice that, in 

the two instances above, the values of Q /Q are quite different. 

This implies radically different turbulence producing mechanisms, 

an implication which is qualitatively consistent with our 

intuition. However, it remains to explore similarly obtained 

2 

values of Q /Q and (Av/u.) for other unstable systems before 
mx o d 

new quantitative insight can be derived. 




H 

co 


1 


Log I 
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IV. Driver Gas Flow with Fluctuations 

A. Introduction 

The conventional Navier Stokes equations do not adequately 

describe turbulent flow. None of the phenomenological approaches 

which have successfully characterized some turbulent systems 

seems derivable from a molecular view of matter. On the other 

14,15,16,17 

hand, a recent series of papers suggests. an alternative 

non-empirical view of turbulence which is based on a modification 
of the usual view of molecular chaos. A set of modified Navier 
Stokes equations is derived as well as a variety of unique 
physical consequences. 

Among these is a characterization of transition to turbu- 
lence in boundary layer flow without the usual restriction to 
linear instabilities. Contact surfaces are ordinarly thought to 

be boundary- layer like insofar as transition to turbulence is 
18 

concerned. Linear approaches seem to have failed in many 

19 

important cases. Therefore, we have explored the potential 
usefulness of contact surface instabilities in driver gas flow 
in a shock tube for new insight on turbulence, paying special 
attention to the possible applicability of the new approach 
from kinetic theory to our observations . 

B. Background: Shock Tube Flow 

For idealized shock tube calculations, infinitely fast 

rupture for the diaphragm material separating the high pressure 

and low pressure sections is usually assumed. (See, for example, 
20 

Thompson on one dimensional unsteady flow. ) This provides the 
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well-known flow regions indicated in the (x,t) diagram in 

Figure 2: a region (1) of un-shocked, driven gas at the original 

driven gas conditions (p^, T^, p^); a region ( 2 ) of shock heated 

driven gas, bounded by the primary shock wave (S) , moving at 

velocity w , on the "downstream" or positive x side and the 
s 

contact surface (c) (separating the driven gas from the driver 

gas), moving at velocity v^i on the "upstream" or negative x side; 

a region .(3) of expansion cooled driver gas; and a region (4), 

of uncooled stationary driver gas at the original gas conditions 

(P 4 , T^, P 4 ) • Between regions (3) and (4), an adiabatic 

expansion is bounded on the upstream side by the head of the 

expansion fan moving upstream (i.e., toward negative x) at the 

local speed of sound (a^) and, on the downstream side, by the 

tail of the expansion fan moving at speed (v^-a^) where v^ is the 

local gas velocity (= v 2 ^* When (v^-a^) is positive, the motion 

of the tail is downstream (i.e., toward positive x) ; cooling of 

the driver gas persists into the shock tube region downstream 

of the plane of diaphragm rupture. In the above, the symbols p, 

T, and p are, as usual, pressure, temperature, and density; the 

subscripts denote the region. Velocities cited are taken from a 

laboratory fixed frame of reference. 

! The contact surface (C in Figure 2, ideally plane and 

18 

stable) is inherently unstable against small transverse 

perturbations. The deviations from one-dimensionality which 

turbulence implies in turbulent detonation waves are usually 

thought to be a result from the shear induced by the curvature 

21 

of the discontinuity. Turbulent contact surfaces also show 
20 

curvature. In addition, nonuniformities associated with the 
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rupture of a diaphragm produce perturbations in the contact 
surface. .Depending on the initial driver gas pressure, p^, 
the local Reynolds numb.er at the contact surface can be made 
high enough so as to be consistent with the possibility of 
eventual turbulent flow. Gas particles are accelerated from 
rest and acquire finally a constant downstream velocity v^ at 
the completion of the adiabatic expansion. Furthermore, in 
the contact surface and wherever else the appropriate conditions 
exist transition to turbulence can take place. 
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C. Apparatus 

Our 12.7 cm diameter shock tube's test section has been 
placed immediately downstream of the diaphragm. High purity 
N 2 gas (or N 2 with various low concentration N 2 0 4 ^ 2N0 2 a< ^“ 
mixtures) expands adiabatically in the 153 cm long driver tube, 
after the rupture of a mylar partition, evolving into a low 
pressure region, the 336 cm long driven tube, also containing 
N 2 gas. The test section is 153 cm long. Pressure histories 
are recorded with wall mounted Kistler 606A gauges at the driver 
section's end plate and at five positions in the test section, 
including one positioned at its center. 

In addition, two 2.5 cm windows at the center of the test 

section 76.2 cm (downstream) from the diaphragm allow an optical 

path through the shock tube. In some shock tube firings, this 

path has allowed us to use blue light absorption (0.435 ym) 

from a 200 Watt mercury arc lamp as a diagnostic for NC> 2 density. 

In other firings, the optical path has been used as one. arm of 

a Michelson-Morley interferometer. In both cases, the basic 

22 

set-up is conventional (see, for example, Bradley for a dis- 
cussion of these techniques). However, in the latter instance, 
simultaneous gas density measurements are made possible along 
two lines equidistant from the diaphragm and a distance one cm 
apart transverse to the flow of gas in the shock tube. The 
amplifications and the d-c supply voltages for the two photo- 
tubes which are used as fringe shift sensors are adjusted so 
that, for laminar unsteady low velocity flow, their signals 
are correlated and cancel when combine^. Figure 3 shows a 
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sketch of the test section and a sample of the interf erpmeter 
set-up data where ordinary room air flow (resulting from the 
air-conditioning in our laboratory) has been used. 

The pressure gauges provide start and stop signals for 
three time interval counters. Pressure histories, NC^ 
absorption histories and the output from the interferometer 
(signals Fl, F2 and F = F1-F2) are recorded on oscilloscopes 
with bandwidth > 30 MHz. A Biomation 8100 A to D converter is 
used in selected instances so that some of our data can be 
subjected to computer analyses; either the pressure history at 
the test section window is digitized along with the simultaneous 
optical data or pressure histories at two different locations 
are digitized. 
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D. Results 
(i) Calibration 

We have tested the use of the equations of one-dimensional 

gasdynamics in our experiments by comparing the measured values 

of the shock strength (p 2 /P 1 ) with the values of p 2 /p^ calculated 

using the 1-D equations and the measured primary shock wave 1 s 

Mach number, M (= w /a.. ) . The range 1.6 < M <2.1 is chosen; 

S S X s 

for a N 2 /N 2 system (N 2 as the driver and driven gas) in our 

shock tube, this range has given us the opportunity of varying 

the Reynolds numbers in the driver gas at fixed M g since p 4 can 

be varied independently of P^/p^ (which determines M g ) . Figure 4 

shows the results from these comparisons. On the average, the 

measured values of P 2 /Pg are higher than calculated by an amount 

10 ± 11%; this effect gives an indication of the limit to which 

we might reasonably employ the 1-D equations in order to predict 

flow conditions in the shock wave processed driven gas. 

The usefulness of the interferometer data as an indicator . 

of contact surface arrival has also been tested. We monitor 

the output of one of the two phototubes on which a single fringe 

has been placed. We observe the sudden change in the fringe 

shifts to a condition of large (and apparently highly erratic) 

variations. We define this event as the arrival of the contact 

surface and measure the time interval which has passed since the 

arrival of the shock wave at the same point; this time interval 

is labeled T . For the same shock tube firing, we use the 
c ,m 3 

measured Mach number M to calculate the time interval which 

s 

would be expected according to the equations of one-dimensional 



Ap/p 


22 ' 



“ 1.6 1.8 2.0 


4 



23 


gasdynamics; this time interval is labeled' T . Figure 5 

C f c 

shows comparisons of T c c and‘T c These data fit a straight 

line through the origin (~1 ± 42 psec) with a slope of one 
(0.99 ±-.08) and a correlation coefficient (r = 0.97) which is 
very close to one. Thus, the measured Mach number seems to be 
an even better predictor of the velocity of the contact surface 
than it is of the pressure ratio across the shock front. 

In addition, we have tested the accuracy of blue light 
absorption as a diagnostic for NC^ density in our set-up. 

For these calibrations, absorption was measured with a static 
fill; there was no diaphragm between the driver and driven 
sections of the shock tube. Lambert-Beer 1 s Law was verified in 
the form: 


log (I/I ) = a + b INCU 
o ^ 

where b = ed, e is the extinction coefficient, d is the- path 

length (= 13.97 cm) and [N0 2 3 is the N0 2 density in mole/liter. 

We found a = -(2.9 ± 11) x 10 ^ and s = 149, ± 8 A/mole-cm. (See 

Figure 6.) This value of e is consistent with previous measure- 

23 

ments as reviewed by Zimet. 
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(ii) Fluctuations 

Significant localized deviations from the mean behavior 
are found in data from the pressure measurements, the inter- 
ferometer and the blue light absorption. These deviations 
are found in the driver gas, upstream from the contact surface. 
The pressure data show large localized variations from the 
average local pressure. These are associated with a simultaneous 
change in the behavior of the interferometer output (F) from 
random variations to a behavior which shows that the density 
changes sensed by Fl and F2 are strongly correlated. When a low 
NC >2 concentration mixture is used t^+^O^^+ZNC^) as the driver 
gas, the pressure deviations are also associated with simul- 
taneous large decreases in the blue light absorption corres- 
ponding to decreases in the NO 2 density. Examples of each of 
these are shown in Figures 7 , 8 and 9 . 

The statistical behavior of the pressure history has been 
studied. We use the correlation function R(x) : 

R(t) = p* (t)p' (t+T ) /p 1 (t) 2 , 

where p' (t) is the difference" between the instantaneous 

pressure and the average pressure at time t. When we compare 

pre- and post-fluctuation pressure data with data obtained during 

the fluctuation, the correlation function appears to be 

structured for data in the fluctuation and unstructured for data 

before and after the fluctuation. A sample of this is shown 

in Figure 10. The behavior of R(x) during the fluctuation is 

24 

consistent with the behavior expected in turbulent flow. 
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Other dynamical features of the fluctuations have been 
measured. At a fixed observation station, the width of the 
deviation in the pressure history twp (that is, its duration 
in the laboratory frame of reference as observed on the 
oscilloscope trace) appears to be relatively insensitive to 
increasing time lag tD behind the contact surface. This lack 
of trend is echoed in the duration of the signature for the 
fluctuation in the interferometer data twi. Figure 11a shows 
this for shock tube firings where simultaneous pressure and 
interferometer data have been obtained. Figure lib shows that 
fluctuations become increasingly less noticeable with time; 
the magnitudes of the ratio of the peak pressure in the 
deviation, p g , to the local mean pressure p L , p B /p L , diminish 
with tD ■*" . The speeds of several fluctuations have been measured 
by tracking them through consecutive pressure stations. When 
the measured velocities are compared with the local velocity 
V 3 calculated using the measure primary shock wave ’ s Mach 
number M g , we find that the fluctuations are moving at approxi- 
mately 90% of the velocity of the local flow (AV/V^ = 0.13 ± .12) . 
These data are given in Figure 12. Furthermore, we have displayed 
in Figure 13 a Reynolds number dependence in the time lag behind 
the contact surface tD at which the fluctuations are detected. 
Shown are those shock tube firings in which two distinguishable 
fluctuations are discernible in the pressure history. Hence- 
forth, the first fluctuation will be called the B 1 fluctuation and 
the second called B 2 . The Reynolds number is based on the local 
flow conditions calculated for the fully expanded driver gas 
(region 3) in units of inverse centimeters according to 
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ReVcm = u 3 /v 3 where v is the kinematic viscosity. Notice that 
the separation between the fluctuations seems to be constant for 
all cases, 470 ± 60 ysec. Taken separately, the time lag tD 

4 * 

decreases with increasing Re /cm for both the first (8^) and 
second (8 3 ) fluctuations. 
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(iii) Boundary Layer Flow Model 

Let ' s model the turbulent, contact surface as a boundary 
layer which is transverse to and symmetric about the axis of 
propagation. That is/ its smallest dimension is at the 
center of the tube; its largest dimensions (small with respect 
to the radius of curvature of the front) are at the shock 
tube's walls. This apparent boundary layer is perpendicular 
to the "normal" boundary layer at the shock tube ' s- walls . The 
local Reynolds number in it for the observed fluctuations can 
be defined using: a local characteristic length, 

H* H a_t 
3 T 

where t T is a characteristic turbulence time; a characteristic 
"free stream" transverse boundary layer velocity, 

U d = U 3 < e2 / 8 > 

where v is the total angle subtended by the .curved contact 
surface; and the values of and computed from the 1-D 
gasdynamical equations. From Figure lib, we also assume that 
the turbulent intensity (5 Q) is related to the time lag of 
the fluctuations as observed at our test section according to 

Q a P B /P L and 

Q = f ± (i,...) (tDB i )“ 1 IV-(l) 

In this, i is a fluctuation ordering parameter and. f is - some 
(unknown) function of i whose overall value is constant in 
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our data once the parameter i is fixed. This assumption does 
not eliminate the possibility of a M s or dependence; in the 
data for which the above relationship is confirmed (i.e., the 
data in Figure 11b)/ the variation in M was only ±10%. For 
consistency, our data require since we find, when 

tnq = toe 2 , (P B /p L ) 6l < <yp L )(s 2 - 

In our case, t is determined from the correlation functions 

in our data, as illustrated in Figure 9, to be t ~ 40 ysec; 

19 24 

0 is estimated from the data reported in Hall and Thompson 
to be 0= 0.07 ± 0.02 rad. The data displayed in Figure 13 
(tD vs Re/cm) is replotted in Figure 14 as log (i^) vs Re where 


Re 


D 3 0 a 3 t T 


In both cases, a good fit to a straight line is obtained. For 

the fit log (l/tD(3 i ) = b^+m^ Re^, we found b^ = -3.8, = -3.7, 

-3 -3 

ny = 3.3 x 10 , and = 2,1 x 10 . Although the correlation 

coefficient r is nearly one in both cases , the large uncertainty 
in 0 quoted above leaves the fit parameters very uncertain 
('v ±50%). Nevertheless, Figure 14 suggests quite clearly a 
linear dependence of the logarithm of turbulent intensity on 
local Reynolds. 

For transitional boundary layer flow, Tsuge's (1974) 
solutions dictate the familiar form 


I E Q/Q q = exp [<j>(Re,A)3 IV- (2) 

The - expansion parameter X , determined by boundary conditions , 
has the units of frequency; Q q is a minimum intensity level which 
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must be exceeded in order for fluctuations to grow. The 

function <j> has no arbitrary constants. The family of curves 

of I* vs Re (where I‘ = log I) at selected values of the 

2 

nondimens ional ratio Av/U^ shows: a first regime with 

(dl'/ciRe) > 0, (d^I'/clRe^) > 0; a middle regime which is 

essentially linear, with (dl’/dRe) > 0, (d 2 /I'/dRe 2 ) = 0; 

and a third regime with (dl'/clRe) > 0 initially and 
2 2 

(d I'/dRe ) <0. The linear middle regime has a slope which 

we shall name m and an intercept on the Re axis which we 

name Re^. Extrapolated to I* = 0, the middle regime gives 

at Re = Re^, Q = Q 0 * T ^e third regime has a maximum in 1 1 of 

value I ' at a Reynolds number which we now name Re . For 
rax p 

—3 

the Reynolds number range in Figure 13, in = 3 x 10 and 
(Re^-Re n )= (0. 56) Re^-0.29. Since, in the linear regime, the 
experimental data should fit log Q = b + (m) (Re) , the deter- 
mination of Re gives an estimate of Q from: 
n J o 

log Q l Re=Re = log Q q = b + (m) (Re n ) IV- (3) 
n 

Returning to Figure 14, and testing the approach above, 
one notices that both of the slopes are consistent with the 
theoretical value within the uncertainty of our estimates . 
There is no evidence of a peak in the experimental Q vs Re 
profile. However, since Figure 14 tells us that Re^ > 310 
if we assume that the theory is applicable, then Re^ > 137 
from the discussion above. Combining this inference with 
equations IV- (1) and IV- (3) and the fit parameters for the 3^ 
fluctuations we find 

(Q/Q q ) > 4. 
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E. Conclusion 

Our shock tube configuration has provided us access to 
gasdynamical phenomena not previously reported. From three 
different diagnostics, viz., pressure measurements, total 
gas density detection, and reactant density detection, we 
find localized fluctuations in the gasdynamic variables in 
the driver gas flow which are distinguishable from the average 
behaviors. These fluctuations show systematic changes with 
Reynolds number as well as unique and reproducible features 
in the nature of their characterizing parameters , viz . , time 
lag with respect to the contact surface, duration, deviation 
strength, and statistical correlation time. These trends do 
not allow the interpretation of our fluctuations as a spurious 
wall boundary layer phenomena. Neither do these trends seem 
affected by the variety of diaphragm materials and diaphragm 
rupturing procedures used. However, the trends are not 
inconsistent with an interpretation of the observed fluctuations 
as examples of turbulent bursts, originating in. an unstable 
contact surface. In addition, attempting to interpret the 
contact surface as a boundary layer in accordance with Tsuge's 
kinetic theory of turbulence, we find a qualitative success 
in terms of the Reynolds number trend with turbulent intensity. 
We can compare the form taken by our data with the form taken 
by theoretical solutions in the same Reynolds number range and 
estimate that the ratio of required maximum to minimum turbulent 
fluctuation level for amplifiable instabilities must be at 
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least four. This overall approach should therefore be extended 
to and tested on a wider class of unstable gasdynamic discontin- 
uities. 
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y. Improving the .Mach Number Capabilities of Arc Driven Shock 
Tubes . 

A. Introduction 

25 

The pressure loaded arc driven shock tube has proven to 
be an invaluable tool in the study of high speed, high 
temperature fluid phenomena and in the creation of high 
temperature plasmas. The continuing need for improved overall 
effectiveness for such devices has motivated a lot of work on 
the effects of circuit inductance, capacitance, impedance 
matching and driver tube configuration on the efficiency of 

25-29 

conversion from electrical energy to test gas kinetic energy. 

The upper limit of Mach numbers attainable from shock tubes has 
correspondingly increased. However, insofar as systematic 
trends in performance estimators are concerned, these previous 
results have seemed to be either tentative or inconclusive. 

• We, therefore, have analyzed the scaling of shock tube 
performance with capacitor energy, initial' driver gas pressure 
and driver volume. We have used data published in the literature 
from ambitious and large configurations as well as the measure- 
ments from our own modest facility. Specifically, we have 
chosen to try to learn as much as possible from the simplest 
model of the shock wave pressurizing process. With this kind of 
device, we expect to ultimately produce the well-known- collisional 
turbulent shock fronts arising from electron recombination processes. 
Section VI gives background for our planned extensions of current 
diagnostic techniques for application to the shock tube environment. 
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B. The Pressurization Process 

In general, pressure loaded arc driven shock tubes can be 
thought of as operating in the same way as ordinary pressure 
driven shock tubes . The electric arc resulting from a spark 
discharge between two electrodes heats up the driver gas, 
uniformly increasing its pressure and temperature at constant 
volume. The increased pressure and driver gas speed of sound 
should routinely allow greater Mach numbers to be 'obtained in 
the driven tube test gas after bursting a diaphragm separating 
the two sections than would otherwise be obtained without 
preheating. We shall assume that this model adequately 
characterizes the pressure loaded arc driven shock tube. 

The performance of the shock tube is often discussed in 
terms of the highest Mach numbers which it might provide. 

However, the speed of the shock wave obviously depends on 
both the driver and driven gas parameters . Therefore , one should 
not expect to successfully use the Mach number alone when comparing 
the performance of different shock tubes if there are varying 
configurations and operating conditions. 

A more useful basis for comparison can be found as follows. 

30 

Recall, for a strong shock wave, that the Mach number is 
expressed as 

_ T l +1 , T 4 T 4>‘l > 1 '' 2 u . 

2 l Y 1 T 1 V 4 ' c 4 V-CL) 

where u/c^ is an implicit solution of the equation - 




V-(2) 
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.Here, using conventional notation, y is the ratio of specific 
heats, T is the temperature, p the molecular weight, p the 
gas density, u -is the gas flow velocity and c is the speed of 
sound. When this formulation is applied to our model, the 
subscript 1 refers to the quiescent driven gas and the subscript 
4 refers to the driver gas after the discharge but before the 
rupture of the diaphragm. 

The solution of Equation V-(2), u/c^ = F(p 4 /p 1 ), can be 
obtained numerically. Figure 15 shows solutions for two sets 
of y's. We note that F(p 4 /p 1 ) is not a very strong function of 
^ 4 ^ 1 ’ The more important parameter for increased Mach number 
is, therefore, T^, which depends only on driver gas conditions. 

It is well known that the cost of a shock tube facility is 
usually a very sensitive function of the cost of capacitors and 

o 

related power supplies, i.e., on E = (1/2) CV, the energy 

c c 

•available to the driver gas. Thus we have chosen, as a result 
of these considerations , to use the ratio of stored capacitor 
energy E c to shock heated driver gas temperature T^ : 

E 

T1 E ■— Joules/°K V~(3) 

X 4 

as our measure for comparison. The g is an index of performance 
loss since an improvement in energy transfer to the driver gas 
shows up as a reduction in n . 
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C. Performance Analyses 

The most interesting configurations for these arc driven 

shock tubes have been either those with cylindrical drivers 

with exploding wire triggers or those with conical inserts in 

29 

the driver tubes. To us, the range of data available seemed 
deficient in the low energy and small driver gas volume regime. 
Our configuration has been chosen so as to begin to -remedy 
this deficiency. 

in our set-up, we have used a cylindrical driver (with 
electrodes similar to that of the conical driver) , 10 cm long 
with a diameter of 5 cm. Arc heating results when a 12 yP 
capacitor charged to 20 kV is discharged through a spark gap 
in the driver tube to which the capacitor is coupled by a 
parallel plate transmission line. The luminous shock front 
is monitored using piezoelectric pressure transducers, photo- 
multipliers and the usual associated diagnostic electronics. 
Helium and argon were used as the driver and driven gases 
respectively . 

(i) Dependence on E and V. 

.The dependence of the shock speed on driver energy density 
has been computed for several facilities. The results are given 
in Figure 16. The speed shown there is in arbitrary units. Each 
set of data has been multiplied by a constant uniquely determined 
for that set; the constant was determined so that all of the 
data would fall near a single line for the purpose of looking 
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for a simple overall secular trend. It is quite evident that 

1/3 

such a trend exists; the speed scales' as (E /V.) , the cube 

C tr 

root of the ratio of stored capacitor energy to driver volume. 
2 

Since M- T^ and T^ = E c /ri, Figure 16 tells us that 


„ „ e 1/3 v/ 3 
c 4 


V-(4). 


This scaling of M with E may be due to increased losses 
at higher currents and to reduced discharge resistance at high 
temperature (current) . The cost per unit temperature increase 
of driver gas increases with E c . Notice that there is no sign 
of a levelling off of speed with increasing E c even at the 
highest energies. Thus, it may be that the speed one obtains 
depends only on one's willingness to invest in bigger capacitor 
banks . 

Figure 16 also suggests decreasing as a way of increasing 
performance without increasing the cost. Of course, this would 
be achieved at the risk of reduced time available for experimental 
measurements at a fixed test station location. 

An unadjusted numerical comparison of some of the over- 
lapping data is given in Table I. At comparable Mach numbers, 
one sees that increased is associated with increased g. 

Our driver, as would be .expected, has approximately the same 
performance as a conical tube if a scaling of our results to 
the higher energies is performed. 
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(ii) Dependence on Initial Driver Gas Pressure p^ 

For a partially ionized plasma discharge column', the 
resistivity will increase with increasing gas pressure. 
-Increasing the resistance of the gap will increase the energy 
dissipated in the arc. A higher pressure, on the other hand, 
dissipated in the arc. A higher pressure, on the other hand, 
will mean a greater mass of driver gas, a larger driver gas 
heat capacity, i.e., a reduced , and therefore an increased 
index of performance loss q. 

We have measured the dependence of Mach number M on 
The results are shown in Figure 17 . Note the existence of a 
plateau at high pressure where M becomes effectively independent 
of Since in our set up then p^ = RT^ = RT-^. 

The fractional efficiency of electrical energy conversion is 
given by 

e = I f>4 E V ( W 

and we find, using equation V-(l), 

M -v (-l_) 1 / 2 p(p 4 /p i ) 
p 4i * J ~ 

Thus, the plateau in Figure 17 implies that the resistivity does 
not increase with pressure fast enough to compensate for the ’ 
increased heat capacity of the driver gas . 

This inference is confirmed in Figure 18. The measured 
efficiency is explicitly shown to have a dependence on p^ which 
is too weak for an increase in p^ to be reliably impactful on 


the Mach number. 
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Thus, it seems unnecessary to' operate the pressure loaded 

arc driven shock tube at the ultra high driver pressures which 

are now often used unless there is some motivation other than 

improved performance. Reduced driver pressure, wherever . 

possible, will eliminate the need for heavy and expensive 

driver tube designs. The lower limit for p 4 ^ should be 

determined, insofar as these analyses are concerned, by the 

31 

need for reproducible flow properties. 

(iii) Dependence on Discharge Length 

The total resistance of a discharge is directly propor- 
tional to its length. For maximum energy dissipation in the 
arc, the length of the arc should be maximized for a given 
driver volume. Using the efficiency as defined above, 

T 4 /T x = 1 + (f £E c /V 4 )/(p 4 RT 1 ) . 

At large Mach numbers, 

t 4 /t x = (2 e c )/(3v 4 p 4 rt 1 ) = E c /nT . 

Hence 

E -1/3 

e - (3V 4 P 4 R/2n) * (—■) V- (5) 

and we see that the efficiency of conversion increases with 

. . 1/3 

increasing discharge length (o, v 4 ' ) . 

Nonetheless, increasing V 4 to increase resistance and 

conversion efficiency is not advisable since it decreases the 

total energy density and increases the index of performance 
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loss T| as discussed above. Attempts, therefore, to improve 
performance through the management of driver gas circuit 
resistance inevitably encounter a fundamental handicap. 

D. Conclusions 

A simple physical model of shock wave pressurization has 
apparently succeeded in providing a consistent overall inter- 
pretation of the performance of a variety of pressure loaded 
arc driven shock tubes. This model should be limited only 

in the sense that the equations of one-dimensional gas dynamics - 

30 

are usually limited in shock tube analyses. Within these 
limits: The change in electrical conversion efficiency resulting 

from a change in the initial driver gas pressure does not 
significantly improve either the performance or the highest 
Mach number available from a given facility once a nominal 
operating p^ has been ascertained. Secondly, increasing the 
discharge volume so as to increase the conversion efficiency 
represents a change in configuration in this model. Table I, 
equation V-(4) and equation V-{5) show that such a change can 
be expected to have a complicated effect on the index of perfor- 
mance loss (and cost) and Mach number; the driver volume should 
be the smallest compatible with experiment test times and the 
discharge length should be maximized for this volume. However, 
at constant volume, Mach numbers and the index of performance* 
loss can be expected to increase (while the efficiency decreases) 
indefinitely with the cube root of capacitor energy. 



TABLE I 


TYPE 

E 

C 

^4i 

V 4 (cm 3 ) 

M 

V T i 

P 1 

■n 

Ref 







air 




cylindrical driver? 
spiral trigger wire 

1 

MJ 

27.2 atm 

5969 

44 

44 

1 

Torr 



1 



He 


34 

41 

10 

Torr 

7 

MJ/K 







22 

37 

100 

Torr 




cylindrical driver; 

1 

MJ 

27.2 atm 

5969 

30 

21 

1 

Torr 



1 

straight trigger 

• 


He 


23 

19 

10 

Torr 

15 

MJ/K 


wire 





16 

20 

100 

Torr 




conical driver 

235 

kJ 

9 . 2 atm 
He 

632 

55 

39 

57 

42 

0.1 

1 

Torr 

Torr 

1.4 

MJ/K 

2 

conical driver 

290 

kJ 

11.8 atm 

350 

85 

132 

0.05 

Torr 

0.7 

MJ/K 

3 




He 


80 

68 

0.05 

Torr* 

1.3 

MJ/K 


Present work 
cylindrical driver 

2.4 

kJ 

less than 

203 

14 

8 

1 

Torr* 




spark gap triggered 



1 atm He 


19 

8 

. 0.1 

Torr* 

0.09 

MJ/K 



* Driven gas is argon. 
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VI. Resonant Absorption in an Argon Plasma at Thermal Equilibrium 

-j ■ 

A. Introduction 

The absorption of resonant argon ion laser radiation has 

been proposed as a sensitive temperature diagnostic for an 

32 

argon plasma in thermal equilibrium. Absorption as a diagnostic 
for electron density and temperature is a well established technique 
in plasma physics that supplements other well known techniques 
such as line and continuum intensity measurements, probe measure- 
ments, light scattering, interferometry, and so on. The main 
advantage of absorption is in its relative simplicity; all that 
is needed is a light source (laser) and a detector without the 
need for calibration. On the other hand, absorption coefficients 
depend on both the density and the temperature and, their use has 
the disadvantage of requiring an independent measurement of one 
of these parameters. The calculation of the absorption coefficient . 
in Ref. 32 has two noteworthy features: (1) the absorptions due 

to free-free transitions (inverse bremsstrahlung) and bound-free 
transitions (photo-ionization) are neglected; (2) the effects 
of Doppler broadening at low densities and high temperatures are 
not included. In this paper, we present calculations of the 

O 

absorption coefficient of an equilibrium argon plasma at X = 4880 A 
without the limitations just mentioned and evaluate their impli- 


cations . 
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B. Theoretical Approach 
(i) The Absorption Coefficient 

The total absorption coefficient K can be expressed as 

K = (1- (v 2 /v 2 ) )~2 (K . + K + K + K ) VI- (1) 

p' ei en p r 

where v is the plasma frequency and v is the laser frequency. 

P » 

The contributions to K .indicated in equation VI- (1) come from 
the following contributions: (i) electron inverse bremsstrahlung 

in the presence of an ion, I (ii) inverse bremsstrahlung by 
electrons during collisions with neutral atoms, K en ; (iii) photo- 
ionization of highly excited atoms and ions, K ; (iv) resonance 

XT 

O 

absorption of 4880 A radiation by Ar II ions, K^. Let's now 

explicitly characterize these contributions. 

The absorption coefficient due to electron inverse 

33 

bremsstrahlung in the field of ions is: 


3.692 x 10 n 

vV '' 2 


(nn + 4n 2 G 2 } 


(c.g.s.) VI~(2) 


where n^, n^ and n^ are the electron, singly charged ion and 

doubly charged ion densities respectively, T is the temperature 

and G. and G„ are the Gaunt factors. The value of G 1 is taken 
x 34 

from Schnapauffs work (Gi = 1.7) and G^ is assumed to be equal 
to one; the Gaunt factors are relatively insensitive to the 
temperature . 

The electron-neutral' atom inverse bremsstrahlung absorption 
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coefficient can be approximated by the expression used in 

35 


microwave theory: 


K 


en 


-(■ 


2 

e n 


7T m c 
e 


v 


ef f 


{v 2 + 


VI- (3) 


where v = n <vcr, > . Here, n is the density of neutrals 
eff o tr o J 

and <va^ r > is the product of the electron velocity and the cross- 

section for momentum transfer collisions averaged over a 

Maxwellian distribution. The value of a. used in the calculation 

36 tr 

is- taken from Engelhardt and Phelps. The validity of using 

equation VI- (3) in the visible ■ wavelengths is discussed in 

37 

numerous papers dealing with laser induced gas breakdown. We 

notice in passing that the electron inverse bremsstrahlung 

absorption during collisions with neutrals is generally much 

smaller than the electron inverse bremsstrahlung except when 

the degree of ionization is less than 1%. 

The contribution of photoionization of highly excited 
33 

atoms is : 


K = 
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2C 1 T , n o*l 
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VI- (4) 
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8.915 x 10 


(c.g.s. ) 
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where g° and g° are the statistical weights of the ground states. 

of the first and second ionization stages respectively and Uq and 

• are the partition functions of Ar I and Ar IX with ionization 

potentials x^_ and x 2 - The lowering of the ionization potentials 

due to the plasma charges are Ax-^ and Ax 2 * The factors are. 

essentially independent of temperature in the visible region of 

38 

the spectrum. The values for obtained by Schliiter and the 

1 39 

values of £ 2 as measured by Vorpaeva are used in our calculations. 

The absorption coefficient of radiation in resonance with 
this transition is 


K r = cP~ B lu ,f,(v) (n H - ^ = IP * (v > n d d-exp (-hv/KT) VI- (5) 

5 U 


where B^ is the Einstein B coefficient and <j) (v) is the absorption 
lineshape. In this, we define n^, n^ and g , g u as the densities 

and statistical weights of singly charged ions in the lower 

2 2 0 
(4s ^ 2 / 2 ^ an< ^ u PP er (4p D 5/2^ levels of the 4880 A argon ion 

transition, i.e.. 


g £ g u 

n £ = ^ exp (-E^/KT) , ^ n ± exp(-E u /KT). VI- (6) 

The reduction in absorption due to stimulated emission is incor- 
porated in equation VI- (5) and the laser linewidth is assumed to 
be much smaller than the absorption linewidth. 

The lineshape factor <f> (v) is a convolution of the (Gaussian or) 
Doppler and Stark broadened lineshapes . Doppler broadening will 
dominate in the low density, high temperature case while Stark 
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broadening will dominate at high densities. A large fraction of 
the total Stark broadening will be due to electron impacts which 
produce a Lorentzian lineshape about a central frequency, v q : 


(Av /2ir) 

4» s (v) - 2 ~ 2 

(v-v o ) 2 + (Av s/2 ) 2 


VI- (7) 


where Av g is the full width at half intensity maximum. ' 

In addition to Stark broadening, there is a Stark shift of 
the line center, so that v and the laser frequency v do not 

O -L* 

coincide, of magnitude d where 


v = v T + d VI- (8) 

o L 

Neglecting ion broadening, Av and d will be proportional to n 

S 6 

and are weak functions of the temperature. (The validity of the 
above assumptions concerning Stark effects and the errors involved 
therein are fully discussed in Ref. 32.) The linewidth and the 

O 

Stark shift parameters for the 4880 A transition are taken from 

40 

the measurements by Labat. 

At very high electron densities, other Ar I and Ar II lines 
may be sufficiently broadened and shifted so as to contribute to 
the absorption. However, since this only occurs when the 
bremsstrahlung and photoionization dominate, resonance absorption 
at the far wings of other lines is not included in our calculations 
given below. 
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(iil Plasma Composition 

For a fixed electron or gas density, the plasma composition 

,is computed using the Saha equation and the laws of conservation 

of mass and charge. We satisfy conditions for local thermo- 

41 

dynamic equilibrium as discussed by Griem. The lowering of 

the ionization potentials is determined from the Debye-Huckel 
35 

theory : 



VI- (9) 


The partition functions are computed by summing over the energy 

42 

levels tabulated by Moore and integrating over assumed hydro- 
genic levels for the excited states close to the continuum. The 
same ionization potential lowering is assumed in the calculation 
of the partition functions and in the Saha equation. Plasma 
composition is computed through the first, second, and third 
ionization stages h.e., for n-^, n^, and n^) • 

(iii) Limitations of Theory 

The limits on the laser - power to be used, W, have been 
discussed in Ref. 32. The lower limit is determined by the ratio 
of the transmitted laser power to the spontaneous emission. The 
upper limit is determined by the condition that the absorption of 
laser energy should not appreciably perturb the system. The 
result in Ref. 32 is: 

2 x 10“ 21 n - W - 4 x 10 -32 n 2 
e e 

-5 

for a beam diameter of 2 mm and a beam divergence of 10 sterad. 

The criterion for the existence of LTE has been formulated 
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by Griem. To insure collisional domination up to the ground 

state, the radiative population rate of the ion ground state from 

the resonance level should be 10 times smaller than the corres- 

32 

ponding collisional population rate. This translates to 


n 

e 


> 


2 x 10 


13 t 1/2 


-3 

cm 


12 

Furthermore, calculations based on a collisional radiative model 

indicate that at n = 10^ cm 3 LTE exists for the Ar I and II 

e 

systems at an electron temperature of T g = 4 eV. Since the lower 

the temperature the easier LTE is reached, then LTE is assured 

when n - 10^ cm 3 and T - 25,000 K. 
e e 

The assumption of Stark broadening by electron collisions 

only entails an error of about 8% for 10^ 6 - n - 10 ^ cm 3 and 
< < 32 

10,000 - T - 40,000 K. Finally, the assumption of a very narrow 
laser line compared to the Stark broadened absorption line is 
easy to satisfy with commercial argon ion lasers especially if 
one uses a single mode laser. 

C. Applications: 10*^ cm' 3 <n <10"*'^ cm 3 ; 10^ °K <T <2.4(10^) K 

Using equation VI- (1), the absorption coefficient K is first 

computed as a function of temperature at fixed gas densities 

n (=n + n, + n„ + n_). The results are shown in Figure 19. Foe 

o 1 2 3 

a fixed n, we see that absorption is sensitive to temperature. 
However, for n ~ 10 17 cm -3 , the absorption coefficient is not a 
single valued function of the temperature. This is due to a 
shift in the balance between the free-free and bound-free 
transitions’ contributions to the absorption coefficient (which 
decrease with increasing temperature) and the contribution from 
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resonance absorption (which increases with increasing temperature) . 

Figure 20 shows our results when the absorption coefficient 
is computed as a function of temperature at fixed electron 
densities. Except at very low electron densities (where the 
validity of an LTE assumption becomes marginal) , the absorption 
coefficient is not very sensitive to temperature in this case. 

In addition, the absorption coefficient is a multi-valued function 

♦ s 

of temperature throughout the computed range. This result is 
markedly different from those in Ref. 32 where no dependence of 
the absorption 'coefficient on electron density was obtained. 

We, therefore, also calculated the dependence of the 
absorption coefficient on the electron density at fixed temperatures 
Here, K is very sensitive to temperature. However, the behavior 
becomes complicated, with increasing temperature, by a plateau 
regime in which a large change in the electron density can leave 
the absorption coefficient relatively unaffected. (See Figure 21.) 
Notice that the dependence of the absorption coefficient on 

electron density is reliably monotonic, in any event, when 

^, A 18 -3 

n >10 cm ■ . 
e 

D. Conclusions 

Of course, the range of temperatures within which accurate 
absorption measurements can be made will depend on the absorption 
length. However, we have found that neglecting free— free and 
bound-free absorption at high electron densities leads to large 
errors in the theoretical estimates of the relationship between^-,., 
absorption and temperature, viz., the serpentine trends shown in 
Figures 19 and 20 are missed. In addition, neglecting the effects 
of Doppler broadening leads incorrectly to the conclusion that 
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the absorption coefficient is insensitive to electron density. 

32 

Therefore, the earlier suggestion that absorption be used to 
obtain temperature measurements rests on an inappropriately 
oversimplified theoretical approach. From our results, it 

O 

seems that absorption of 4880 A radiation by an argon plasma 
(once the temperature is determined independently) is more 
suited as an electron density diagnostic. 
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VII. Summary 

Tests for the existence and dynamic features of turbulent 
collisional shock fronts have not yet been performed. However, 
systematic trends in the performance parameters of pressure 
loaded arc driven shock tubes have been determined. The initial 
driver gas pressure is relatively unimportant. Adjusting the 
driver colume and the length of the discharge path can produce 
mixed results. Nonetheless, under the specified conditions, 
the Mach number increases with the cube root of capacitor 
energy without apparent limit. The relevance of fluorescence 
to ion density measurements in the ionizing shock wave has been 
confirmed. The absorption coefficient of an equilibrium argon 
plasma is computed, avoiding earlier restrictive assumptions, 
as a function of temperature and electron density. Our results 
suggest that the diagnostic applications of resonant absorption 
to temperature estimates are more difficult than previously 
expected. 

When Tsuge's kinetic theory of turbulence is applied to 
flow with a simple nonequilibrium process, we find that the 
reaction rates for the turbulent and nonturbulent cases will be 
apparently different in direct proportion to the degree of 
chaos violation. Our studies of 2N0 2 + N 2 N 2°4 + N 2 an 
adiabatically expanding shock tube 1 s driver gas have revealed 
a systematic Reynolds number based association between 
pressure anamolieS' and distortions in the N0 2 density history. 

We have analyzed these results using Tsuge's theory with its 
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boundary layer approximation. This approach seems to provide 
us with a nonempirical explanation of the well-known "hot spots" 
in turbulent nonequilibrium systems. 

A boundary layer model allows us to determine characteristic 
scales and Reynolds numbers for a variety of turbulent shock 
tube discontinuities. With this model, Tsuge's kinetic theory 
for boundary layer flow with transition to turbulence becomes 
applicable. Minimum fluctuation levels as well as characteristic 
frequencies can be determined for unstable boundary layers 
with turbulent bursts . Previously published results on detonation 
waves have been successfully treated with these techniques. 

In addition, we find that our shock tube's driver gas flow 
produces fluctuations which can be successfully interpreted as 
turbulent bursts in contact surface flow. Therefore, our 
data and analyses support the approach to the theory of turbulence 
in Tsuge's works 'and, in particular, suggest a possible molecular 
interpretation for large coherent structures. 
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